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Abstract 

Basic ingredients of the theory of baryogenesis are reviewed with emphasis on 
out-of-equihbrium decays of heavy particles. The present use of kinetic theory 
is explained and some attempts to go beyond the classical Boltzmann equations 
are discussed. 
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1 Scenarios for baryogenesis 



One of the main successes of the standard early-universe cosmology is the prediction 
of the abundances of the light elements, D, ^He, '^He and ^Li. Agreement between 
theory and observation is obtained for a certain range of the parameter t], the ratio 
of baryon density and photon density , 



where the present number density of photons is ~ 400/cm'^. Since no significant 
amount of antimatter is observed in the universe, the baryon density yields directly 
the cosmological baryon asymmetry, ub — ris — 

A matter-antimatter asymmetry can be dynamically generated in an expanding 
universe if the particle interactions and the cosmological evolution satisfy Sakharov's 
conditions 0, i.e. 

• baryon number violation 

• C and CP violation 

• deviation from thermal equilibrium . 

Although the baryon asymmetry is just a single number, it provides an important 
relationship between the standard model of cosmology, i.e. the expanding universe 
with Robertson- Walker metric, and the standard model of particle physics as well as 
its extensions. 

At present there are a number of viable scenarios for baryogenesis. They can be 
classified according to the different ways in which Sakharov's conditions are realized. 
Already in the standard model C and CP are not conserved. Also baryon number 
(B) and lepton number (L) are violated by instanton processes In grand unified 
theories B and L are broken by the interactions of gauge bosons and leptoquarks. 
This is the basis of the classical GUT baryogenesis Analogously, the L violat- 
ing decays of heavy Majorana neutrinos lead to leptogenesis 0]. In supersymmetric 
theories the existence of approximately flat directions in the scalar potential leads 
to new possibilities. Coherent oscillations of scalar fields may then generate large 
asymmetries f^. 

The crucial departure from thermal equilibrium can also be realized in several 
ways. One possibility is a sufficiently strong first-order electroweak phase transition 
[|^. In this case CP violating interactions of the standard model or its supersymmetric 
extension could in principle generate the observed baryon asymmetry. However, due 
to the rather large lower bound on the Higgs boson mass of about 105 GeV, which 
is imposed by the LEP experiments, this interesting possibility is now restricted to a 
very small range of parameters in the supersymmetric standard model. In the case 
of the Affieck-Dine scenario the baryon asymmetry is generated at the end of an 
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inflationary period as a coherent effect of scalar fields which leads to an asymmetry 
between quarks and antiquarks after reheating [§]. For the classical GUT baryogenesis 
and for leptogenesis the departure from thermal equilibrium is due to the deviation 
of the number density of the decaying heavy particles from the equilibrium number 
density. How strong this deviation from thermal equilibrium is depends on the lifetime 
of the decaying heavy particles and the cosmological evolution. Further scenarios for 
baryogenesis are described in 

The theory of baryogenesis involves non-perturbative aspects of quantum field 
theory and also non-equilibrium statistical field theory, in particular the theory of 
phase transitions and kinetic theory. A crucial ingredient is the connection between 




Figure 1: One of the 12-fermion processes which are in thermal equilibrium in the 
high-temperature phase of the standard model. 

baryon number and lepton number in the high-temperature, symmetric phase of the 
standard model. Due to the chiral nature of the weak interactions B and L are not 
conserved. At zero temperature this has no observable effect due to the smallness of 
the weak coupling. However, as the temperature approaches the critical temperature 
Tew of the electroweak transition, B and L violating processes come into thermal 
equilibrium |]10 . 



The rate of these processes is related to the free energy of sphaleron-type field 
configurations which carry topological charge. In the standard model they lead to an 
effective interaction of all left-handed fermions P] (cf. fig. [^), 

Ob+l = n (iLigLigLiki) , (2) 

i 

which violates baryon and lepton number by three units, 

AB = AL = 3. (3) 
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The evaluation of the sphaleron rate in the symmetric high temperature phase is a 
challenging problem |jTl|. Although a complete theoretical understanding has not yet 



been achieved, it is generally believed that B and L violating processes are in thermal 
equilibrium for temperatures in the range 

Tew ~ 100 GeV <T < Tsph ~ 10^^ GeV . (4) 

Sphaleron processes have a profound effect on the generation of the cosmological 
baryon asymmetry, in particular in connection with lepton number violating interac- 
tions between lepton and Higgs fields, 

^AL=2 = ]^fij lli<P C lLj<P + h.C. . (5) 

Such an interaction arises in particular from the exchange of heavy Majorana neu- 
trinos (cf. fig. In the Higgs phase of the standard model, where the Higgs field 
acquires a vacuum expectation value, it gives rise to Majorana masses of the light 
neutrinos Ve-, i^fi and Ur- 




Figure 2: Effective lepton number violating interaction. 

Eq. d^) suggests that any B+L asymmetry generated before the electroweak phase 
transition, i.e., at temperatures T > Tewi will be washed out. However, since only 
left-handed fields couple to sphalerons, a non-zero value of -B + L can persist in the 
high-temperature, symmetric phase if there exists a non- vanishing B — L asymmetry. 
An analysis of the chemical potentials of all particle species in the high-temperature 
phase yields a relation between the baryon asymmetry Yb = {riB — Ub)/ s, where s is 
the entropy density, and the corresponding B — L asymmetry Yb-l, respectively jl^ . 



Yb = C Yb-l = >l ■ (6) 



The number C depends on the other processes which are in thermal equilibrium |T^ 



If these are all standard model interactions one has C = {8N + 4)/ {22N + 13) for N 
generations. If instead of the Yukawa interactions of the right-handed electron the 
AL = 2 interactions @ are in equilibrium one finds C = —2N/ {2N + 3). 



4 



The interplay between the sphaleron processes (fig. ^ and the lepton number 
changing processes (fig. ^ leads to an intriguing relation between neutrino properties 
and the cosmological baryon asymmetry. The decay of heavy Majorana neutrinos can 
quantitatively account for the observed asymmetry. 



2 Heavy particle decays in a thermal bath 

Let us now consider the simplest possibility for a departure from thermal equilibrium, 
the decay of heavy, weakly interacting particles in a thermal bath. To be specific, we 
choose the heavy particle to be a Majorana neutrino = A^*^ which can decay into a 
lepton Higgs pair and also into the CP conjugate state /0 

N N (7) 



In the case of CP violating couplings a lepton asymmetry can be generated in the 
decays of the heavy neutrinos which is then partially transformed into a baryon 
asymmetry |^ by sphaleron processes [|l^ . Compared to other scenarios of baryoge- 
nesis this leptogenesis mechanism has the advantage that, at least in principle, the 
resulting baryon asymmetry is entirely determined by neutrino properties. 




Figure 3: AL = 1 processes: decays and inverse decays of a heavy Majorana neutrino. 

The generation of a baryon asymmetry is an out-of-equilibrium process which is 
generally treated by means of Boltzmann equations. A detailed discussion of the basic 
ideas and some of the subtleties has been given in [Q. The main processes in the ther- 
mal bath are the decays and the inverse decays of the heavy neutrinos (cf. fig. §), and 
the lepton number conserving (AL = 0) and violating (AL = 2) processes (cf. fig. H). 



5 



Figure 4: AL = and AL = 2 lepton Higgs processes. 



In addition there are other processes, in particular those involving the t-quark, which 
are also important in a quantitative analysis [T^JT^. A lepton asymmetry can be 
dynamically generated in an expanding universe if the partial decay widths of the 
heavy neutrino do not respect CP symmetry, 

r(iv-./0) = i(i + e)r, r(iv^[0) = i(i-e)r, (8) 

where F is the total decay width and the parameter e ^ 1 measures the amount of 
CP violation. 

The Boltzmann equations for the number densities of heavy neutrinos {un), lep- 
tons {ni) and antileptons {ni) corresponding to the processes in figs. ^ and ^ are given 
by 

-.AT) 

- N) , (9) 
N) 

- U) , (10) 

N) 

-l^), (11) 
with the reaction rates 

7(iV^/0) = J d^u3fN{pi)\M{N ^l<f))W... (12) 

= j d^^2Mh{p^)U{p^)\M\l<P-^m\... (13) 

Here H is the Hubble parameter, denotes the phase space integration over 

particles in initial and final states. 



an 



N 



dt 

drii 



+ 3Hn 



N 



+ 3Hni 



-7 (AT /0) + 7(^0 

7(A^ 1(f)) - ^ 

+7(r0^/0)-7(/0 
7(A^ ^ [0) - 7(10 ^ 

+7(/0^/0)-7(^0 




Figure 5: Time evolution of the number density to entropy density ratio. At T ^ M 
the system gets out of equilibrium and an asymmetry is produced. 



and 

fi{p) = exp {-pEi{p)) , ni{p) =gi j^T^fiiP) > (15) 



. (27r)3- 

are Boltzmann distribution and number density of particle i = N,l,(j) at temperature 
T = respectively. Ai and Ai' denote the scattering matrix elements of the 

indicated processes at zero temperature; the prime indicates that for the 2 — > 2 
processes the contribution of the intermediate resonance state has been subtracted. 



For simplicity we have used in eqs. (|T2[) and (0) Boltzmann distributions rather 



than Bose-Einstein and Fermi-Dirac distributions, and we have also neglected the 
distribution functions in the final state which is a good approximation for small 
number densities. Subtracting ([Tl|) from ( p!o|) yields the Boltzmann equation for the 
asymmetry ni — nj. 

A typical solution of the Boltzmann equations (^) - ([TT| ) is shown in fig. |^. Here 
the ratios of number densities and entropy density, 

Y. = ^, (16) 
s 

are plotted, which remain constant in an expanding universe in thermal equilibrium. 
A heavy neutrino, which is weakly coupled to the thermal bath, falls out of thermal 
equilibrium at temperatures T ~ M since its decay is too slow to follow the rapidly 
decreasing equilibrium distribution f^ ~ exp(— /3M). This leads to an excess of the 
number density, un > 'ti'n- C'P violating partial decay widths then yield a lepton 
asymmetry which, by means of sphaleron processes, is partially transformed into a 
baryon asymmetry. 

The Boltzmann equations are classical equations for the time evolution of number 
densities. The collision terms, however, are S'-matrix elements which involve quantum 
mechanical interferences of different amplitudes in a crucial manner. Since these 
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scattering matrix elements are evaluated at zero temperature, one may worry to 
what extent the quantum mechanical interferences are affected by interactions with 
the thermal bath. Another subtlety is the separation of the 2 — 2 matrix elements 



into a resonance contribution and remainder , 

\M{l<P ^ = \M'{l<P + \Mres{l<P W , (17) 

where the resonance contribution has the form 

Mres{l(j) ^ l^) oc M{l(j) N)M{N l^y = \M{l(j) ^ A^)|2 . (18) 

The entire effect of baryon number generation crucially depends on this separation. 
The particles which participate in the 2 — > 2 processes are massless, hence their 
distribution functions always coincide with the equilibrium distribution. Only the 
resonances, treated as on-shell particles, fall out of thermal equilibrium and can then 
generate an asymmetry in their decays. General theoretical arguments require cancel- 
lations between these two types of contributions which we illustrate in the following 
with two examples. 

Cancellations in thermal equilibrium 

If all processes, including those which violate baryon number, are in thermal 
equilibrium the baryon asymmetry vanishes. This is a direct consequence of the CPT 
invariance of the theory. 



{B) = Tt{pB) =TT[{CPT){CPT)-^exp{-f3H)B'^ 

= Tr (exp {-(3H){CPTy^B{CPT)) = -Tr(pfi) = . (19) 

Hence, no asymmetry can be generated in equilibrium, and the transition rate which 
determines the change of the asymmetry has to vanish, 

^^^^^^ + 3if(n, - n,) = A7- = , (20) 

where the superscript eq denotes rates evaluated with equilibrium distributions. 
From eqs. (§), (|10D and (pH]) one obtains for the resonance contribution, i.e. decay 



and inverse decay, 

A^^l = -2eY'iN ^ /0) . (21) 

This means in particular that the asymmetry generated in the decay is not compen- 
sated by the effect of inverse decays. On the contrary, both processes contribute the 
same amount. 

The rate Ajres to be compensated by the contribution from 2^2 processes 
which is given by 

A72I2 = 2 / t/$i234/r(Pi)/r(P2) {\M'ii<p - i0)r - \M'm . (22) 
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For weakly coupled heavy neutrinos, i.e. F oc A^M with ^ 1, this compensation 
can be easily shown using the unitarity of the 5'-matrix. 
The sum over states in the unitarity relation, 

E {\Mil<P - X)|2 - \MiX = , (23) 

X 

can be restricted to two-particle states to leading order in the case of weak coupling 
A. This implies for the considered 2^2 processes, 

y:' {\Mii<p ^ m\' - \Mii'^' = , (24) 

where the summation J2' includes momentum integrations under the constraint of 
fixed total momentum. From eqs. ([T7|) and (0) one obtains 

A72I2 = 2 / C?$1234/r(pi)/;'(p2) {-\Mresil<P ^ T*^)!' + IMresiU ^ /</>)!') • (25) 

In the narrow width approximation, i.e. to leading order in A^, this yields the wanted 
result, 

A72I2 = 2 J (i$i234/r(pi)/r(P2) {-\Mii<p ^ N)\'\MiN 

+ \M{U - N)\'\M{N ^ /0)p) ^6{s - M^) 
= 2e7^«(iV ^ /0) = -A7;L . (26) 

This cancellation also illustrates that the Boltzmann equations treat resonances as 
on-shell real particles. Off-shell effects require a different formalism which will be 
discussed in Sec. 3. 



Cancellations at zero temperature 



The lepton asymmetry which is obtained by solving the Boltzmann equations 
depends crucially on the separation of the 2^2 scattering amplitudes into resonance 
contribution and remainder. How to perform this separation appears obvious for a 
graph like fig. |^ which represents the interference between the tree level and a one- 
loop correction term for the vertex AiV/0, together with a s-channel iV-propagator. 



The identification of the resonance contribution shown in fig. |^ is less obvious. 
One expects mixing effects in the case of several heavy neutrinos which have been 
discussed in ||T7| , p!8|Jl9[| . Treating the one- loop self-energy like the one-loop vertex 
correction yields indeed a finite contribution to the CP asymmetry [T3JT^ which is of 



the same order as the vertex contribution. However, one may worry about the fact 



9 




Figure 6: One-loop vertex correction to lepton Higgs scattering. 




Figure 7: One-loop self-energy correction to lepton Higgs scattering. 



that the same procedure yields an infinite result for the total and partial decay widths. 
The self-energy corrections have to be resummed in order to determine the mass of 
the heavy neutrino, its partial decay widths and, in particular, the CP asymmetry. 
It is well known that the properties of unstable particles are defined by position 



and residue of the corresponding poles of the scattering matrix [20|. These poles 
correspond to the poles of the full propagator. 

It turns out that the resummed propagator does not contribute to the CP asym- 
metry of 2 ^ 2 processes at fixed external momenta |^ (fig. |^). Furthermore, even 



N 

o 



/ 2 



N 

o 



/ 2 



Figure 8: Propagator contribution to CP asymmetry for fixed external momenta. 

the total CP asymmetry vanishes to leading order in when integrated over phase 
space 1^ (fig. 1^). This can be seen explicitly in ordinary perturbation theory 
for center-of-mass energies below the resonance region, s <C M^, as well as in the 
resonance region, s ~ M^, which can only be studied after resummation |pT |. 
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Figure 9: Total CP asymmetry to leading non-trivial order in the coupling. 



The cancellation of the various contributions to the CP asymmetry follows directly 
from unitarity. In fact, this is the physical meaning of eq. (0), 

y:' (i-m(/0 - I'n' - mi'^' ^ m^) = o . 

i,4>,i',<t>' 

Away from resonance poles, where ordinary perturbation theory holds, the CP asym- 
metry vanishes to order A®. Corrections due to four-particle intermediate states are 
0{X^). In the resonance region the CP asymmetry vanishes to order with correc- 
tions 0{X^). 

These cancellations between various contributions to the CP asymmetry demon- 
strate the importance of identifying correctly the resonance contribution. This is 
complicated by the fact that different chiral projections of the resummed propagator 



are diagonalized by different unitary matrices |21]. These matrices, together with 



the vertex corrections, determine the effective Nlcj) vertex (fig. TU), where N now 
corresponds to a pole of the full heavy neutrino propagator. 




Figure 10: Effective Nlcj) vertex including mixing effects. 

Given the effective Nlcf) vertex it is straightforward to determine the CP asym- 
metry in the decay of a heavy Majorana neutrino. 



T{Ni ^ let)) - T{Ni ^ 
T{N,~^lcl)) + T{N,^U) 



(27) 



To leading order in A^, the asymmetry is a sum of two terms, a mixing contribution 
(ir = AtA), 



(28) 
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which is directly related to the one- loop self energy |T^JT^, and the familiar vertex 
contribution 

where 

/(x) = v^(^l-(l+x)ln(^^)) . (30) 

These results hold for sufficiently large mass splittings, i.e. |Mj — Mj\ ^ \Ti — Tj\. 
For small mass differences one expects an enhancement of the mixing contribution 



23| . At present, however, the influence of the thermal bath on this enhancement is 
unclear. 

The use of classical Boltzmann equations with collision terms given by S'-matrix 
elements which crucially involve quantum interferences is unsatisfactory. Like the 
collision terms also the time evolution of the system should be treated quantum 
mechanically. 



3 Quantum mechanics of baryogenesis 

One would like to have a full quantum mechanical treatment of baryogenesis. Starting 
with a density matrix with no initial asymmetry the emergence of an asymmetry 
should be seen from the full time evolution. The Boltzmann equations should then 
follow in some limit first-order approximation. 

The Boltzmann equations are an on-shell approximation: between the interaction 
processes the particles propagate on-shell, and for the scattering processes the on-shell 
S'-matrix elements are used. A full quantum treatment will also take off-shell effects 
into account. 

Of course, it is hard to obtain a full quantum description of the processes as 
the whole system is out of equilibrium. We shall therefore discuss a toy model for 
a relaxation process, following Joichi, Matsumoto and Yoshimura |]2^, and compare 
the exact description with the Boltzmann approach. 

The model consists of a single quantum mechanical oscillator coupled to a thermal 
bath of oscillators. The evolution of the system is governed by the Hamiltonian 

/•oo r 

H = Ec^c+ du dn ub^{n,u)b{n,uj) (31) 

J UJc ^ 

+ jliujjdVL ((:{VL,uj)h\VL,uj)c+ h.c.) . (32) 

E is the frequency of the single oscillator, the frequencies of the bath oscillators uj are 
bounded below by ojc with E > Uc- fl is an additional discrete or continuous label 
for the bath oscillators, and ({Q,uj) denotes the complex coupling. 
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E 



Figure 11: One single oscillator with frequency E is coupled to a large thermal hath 
of oscillators. 



The operators c and h fulfil canonical commutation relations, 

[c,ct] = l , [c,W{yi,u)] = Q , (33) 
[6(0, uj), b^{n', ou')] = (5(0 - fl')S{u} - ou') . (34) 

The Hamiltonian is quadratic in c and b{fl,uj) and hence solvable by a change of 
variables. It is possible to give an explicit formula for new operators B{D,,u!) which 
diagonalize the Hamiltonian, 

/•oo r 

cLu dfl ujB\fl,Lu)B{n,uj) . (35) 

These operators B are obtained as a linear combination of c and h. For uj away from 
E it takes the form 

B{^l,u)^h{^l,uj) + 0{Cc,CCh) . (36) 
The time evolution in the new variables is just the free one, 

B\VL,uj,t)=e''''B\VL,uj) . (37) 

By inverting the change of variables one obtains explicit formulas for the time evo- 
lution of c and h. Hence we are able to discuss exactly the properties of the system. 
Let us start with a simple example. 

Decay process 

Assume that the system is prepared in the initial state \ip) = c^|0) where only the 
single oscillator is excited. We expect that this excited state decays. This should be 
seen from the time evolution of the occupation number of the single oscillator which 
is 

W^mtm = \g{t)\' , (38) 
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where 

r 1 

git) = / dujcx(u)- — — --TTTe*'^* . (39) 

Here (j{uo) is the absolute value squared of the coupling summed over the internal 
label VL, 

a{u}) = fdn\C{n,uj)\'^ . (40) 



We see that g is essentially the Fourier transform of the c-propagator including the 
"self energy" n(a;) + i7ia{u). For weak coupling, i.e. n(_E') <^ E and cr{E) E — Uc, 
the integrand in the expression for g has a sharp Breit-Wigner resonance at u = E. 
The contribution from this resonance is 

go{t) oc e-r*/2+^^* , (41) 

where F = 27ca{E). This is the result we would expect from the Boltzmann equations: 
exponential decay with decay rate F. But this is not the only contribution. There is 
a second contribution from the threshold which at large times only decreases with a 
power law, 

The constants k, a parameterize the threshold behaviour for u close to Uc, 

a{uj) = k{uj - tuj" . (43) 

At large times this contribution dominates the resonance contribution. So we are 
led to the interesting result that the exponential Boltzmann decay is only relevant at 
intermediate times, whereas the asymptotic behaviour is described by a power law. 
This behaviour holds generally in field theory p5| . 

Thermal equilibrium 

We will now study another interesting case, the situation of thermal equilibrium. 
The density matrix is then given by 

P = e-''' , /3=^, (44) 
which leads to the usual Bose- Einstein distribution in the variables B, 

{B\n,u;)B(Q',Lj')) = 6(Q - Q')6{uj - u')--^ — . (45) 

By expressing c in terms of the operators B we obtain the expectation value of the 
occupation number of the single oscillator, 

(etc) = Jduaiu) n(i))2 + (,,(^))2 1 • (46) 
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For weak coupling we have again a sharp resonance near uj = E. The contribution 
from this resonance to the occupation number is 

(^'^) - • (47) 

This is the contribution we would expect for a free oscillator. However, at small 
temperatures the resonance contribution is exponentially suppressed and again the 
threshold behaviour becomes more important, yielding 

C(a + l)r(a + l) 



(etc) 



{E - UJ,) 



\2 



-T°+^ , for Uc<^T<^E 



/t .^i"^^A e-^^-r°+^ , for T<^uJr<E. 



(E - iu.y 

Again we obtain a surprizing result: at small temperatures the suppression of the 
occupation number is much weaker than what is expected from the Bose-Einstein 
distribution. 

What are the implications of this result? In it has been argued that the 
power behaviour significantly affects the WIMP abundance. However, the quantita- 
tive importance of this effect requires further investigations p7| , P^ . 

What are the implications for baryogenesis and leptogenesis? How large are the 
errors for the results obtained with standard Boltzmann equations? To answer this 
question we try to apply our simple model to the case of heavy particle decay. For 
simplicity we consider only scalar particles. Let X be a heavy particle that may decay 
into light scalar particles a and b. The idea is to identify the heavy particle X with 
the single oscillator and the decay products a and b with the thermal oscillator bath. 
The interaction is given by a Yukawa like coupling between the three scalar fields, 

^^int = A /(i3x0x0a06 • 

This interaction part should have the same structure as the interaction part of the 
oscillator model. To see how this identification works we expand the fields in Fourier 
modes. As in the oscillator model we enforce that total particle number is conserved, 
i.e. we only take into account those parts of the Hamiltonian which describe decay 
and inverse decay. 



J (27r)3(2cu,)i/2 (27r)3(2cu,)i/2 {27T)%2uj,y/^ ^ ' 
[c\q)ba{kaMkb)+h.c.) 

= J-0^Jdu;dnan,LU,q){c^{q)bin,u;,q)+h.c.) . (48) 
One then reads off the definition of composite operators, 

1/2 
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which are needed for the identification. Here is a label that together with the 

— * — * 

energy cu and the total momentum q determines the momenta ka, k^, e. g. it can be 
chosen as two angles describing the direction of ka- 

These composite operators fulfil commutation relations that differ from the canon- 
ical ones. To apply the oscillator model we try to approximate the commutation 
relation and the kinetic term of the thermal bath as 



[b{fl,uj,q),b^{Q',u;',q')] = {27rf6{q - q')S{uj - uj')S{Q - Q') , 



bath 



d^q 
(2^ 



dujdilu)b\fl,u),q)b{il,(jj,q) . 



(50) 



(51) 



With these adjustments one can compute the decay of the heavy particle X. The 
approximations made above do not affect the behaviour of X substantially and corre- 
spond to a low density approximation. Again we obtain the result of an exponential 
decay followed by a power law behaviour at large times. 

When describing properties of the thermal bath hke the time evolution of an 
asymmetry things become more involved and it is not possible to describe the long- 
term behaviour in the model. But nevertheless it appears possible to compute the 
time derivative of the asymmetry correctly. 

An important question is how to implement CP violation in the oscillator model 
without breaking CPT invariance. This can only be done in an extended version of 
the oscillator model which involves an extra interaction term. In connection with 
leptogenesis this corresponds to effects of the heavy neutrinos N2 and A^3. 
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FINAL ASYMMETRY 
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Figure 12: The final asymmetry as function of r] = ^ for two values of 
0.0001 (sharp resonance) and 0.1 (broad resonance). 



The production rate of the asymmetry can now be calculated. The resonance 
contribution coincides exactly with the rate given in section 2 including both contri- 
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butions (^) and (pOf ) to the CP asymmetry. A detailed analysis of these results is 
given in p9| . 

Despite the various remaining problems it is worthwhile to estimate the size of 
the corrections for the baryogenesis scenario. In such an estimate is given under 
certain assumptions and approximations. The result is shown in fig. |T2|. If the decay 
rate F is of order the Hubble rate H or smaller there is no significant deviation from 
the result obtained with the Boltzmann equations (on-shell result). For T/H > 1, 
where the asymmetry is smaller because the system is closer to equilibrium, the 
asymmetry production is enhanced by off-shell effects. 



4 Conclusions 

At present there still exists a variety of mechanisms which, at least in principle, can 
account for the cosmological baryon asymmetry. Particularly successful is the lepto- 
genesis scenario. Given the experimental indications for neutrino masses it naturally 
explains the observed order of magnitude of the cosmological baryon asymmetry with- 
out any fine tuning of parameters. 

One is therefore led to examine the theoretical basis for current estimates of the 
baryon asymmetry. Although our understanding of the high-temperature symmetric 
phase of the standard model has significantly improved in recent years the quantitative 
description of an out-of-equilibrium process remains a difficult problem. 

Particularly subtle is the use of the classical Boltzmann equations together with 
collision terms derived from ^-matrix elements which involve quantum interferences 
in a crucial manner. A full quantum mechanical treatment which includes the time 
evolution of the system is highly desirable. At present several interesting ideas are 
pursued by different groups but a fully satisfactory solution of the problem still re- 
mains to be found. 
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